Rings sliding on a honeycomb network: Adsorption contours, interactions, and 

assembly of benzene on Cu(lll) 



K. Berland 1 '*, T.L. Einstein 2 't, and P. Hyldgaard^f] 

x '* Department of Microtechnology and Nanoscience, MC2, Chalmers University of Technology, SE-41296 Goteborg, Sweden 
2 ^ Department of Physics, University of Maryland, College Park, Maryland 20742-4111 

(Dated: October 14, 2009) 

Using a van der Waals density functional (vdW-DF) [Phys. Rev. Lett. 92, 246401 (2004)], we 
perform ab initio calculations for the adsorption energy of benzene (Bz) on Cu(lll) as a func- 
tion of lateral position and height. We find that the vdW-DF inclusion of nonlocal correlations 
(responsible for dispersive interactions) changes the relative stability of eight binding-position op- 
tions and increases the binding energy by over an order of magnitude, achieving good agreement 
with experiment. The admolecules can move almost freely along a honeycomb web of "corridors" 
passing between fee and hep hollow sites via bridge sites. Our diffusion barriers (for dilute and 
two condensed adsorbate phases) are consistent with experimental observations. Further vdW-DF 
calculations suggest that the more compact (hexagonal) Bz-overlayer phase, with lattice constant 
a = 6.74 A, is due to direct Bz-Bz vdW attraction, which extends to ~ 8 A. We attribute the 
second, sparser hexagonal Bz phase, with a — 10.24 A, to indirect electronic interactions mediated 
by the metallic surface state on Cu(lll). To support this claim, we use a formal Harris-functional 
approach to evaluate nonperturbationally the asymptotic form of this indirect interaction. Thus, 
we can account well for benzene self-organization on Cu(lll). 

PACS numbers: 68.43.-h,73.20.-r,71.15.Mb,73.90.+f 



I. INTRODUCTION 

With adsorption of atoms and of simple molecules on 
metals relatively well understood, interest has turned to 
the adsorption of organic molecules There are many 
technological motivations for studying adsorption of elec- 
tronically functional organic overlayers, e.g., organic 
light-emitting diodes (OLEDs) offer great promise for 
displays) 2 ^ Acenes and acene derivatives are of particular 
interest, not just because of their importance in organic 
semiconductors but also because of the self-organized 
structures that they form on surfaces. 

Arguably the most dramatic example of the latter 
is the giant regular honeycomb network formed by an- 
thraquinone (AQ) on Cu(lll)PThe pore diameter is un- 
precedentedly large, over 5nm, and each cell encloses over 
200 uncovered Cu surface atoms. Our understanding is 
that the adsorption process is primarily due to van der 
Waals (vdW) interactions, the short-range arrangement 
is due to hydrogen bonding from an aromatic hydrogen 
on one AQ to an oxygen of a carbonyl group on its neigh- 
bor; we believe that the long-range regularity can be at- 
tributed to indirect interactions mediated by the metallic 
(Shockley) surface state centered around the zone center 
of noble metal (111) surfaces, 4 in a mechan ism analo gous 
to that documented for atomic adsorbates! 5 * 6 * 7 * 8 ^ ^ 

In a first attempt to understand the details of the ad- 
sorption and self-organization of acene derivatives, we 
here consider the simplest member of this family, the 
monocyclic aromatic ring benzene (Bz). Adsorption of 
Bz on Cu(lll) displays rem arkable molecula r-adsorption 
dynamics with observation J 11 * 12 * 13 * 14 * 15 * 1 ^ 7 ^ of high mo- 
bility of individual Bz molecules at low temperatures, as 
well as a possibility for the Bz adsorbates to condense into 



two coexisting adsorbate-overlayer phases, CI and C2, 
characterized by hexagonal nets with lattice constants 
10.24 A and 6.74 A, respectively. For Bz, as for AQ, the 
adsorption energy arises primarily from van der Waals 
(vdW) interactions with conduction electrons on the Cu 
surface, with very small orbital changes and substrate- 
adsorbate charge transfers.^ 

The case of Bz on Cu(lll), Fig. 1, provides a forum 
at which to show the power o f a rece ntly-developed vdW 
density functional (vdW-DF) 19 | 2Q | 21 1 to describe molecu- 
lar physisorption on a system with extremely low cor- 
rugation. Systems with low corrugation are interesting 
and challenging since they leave the adsorbates free to 
respond to weak direct and indirect mutual interactions 
which must therefore be described and analyzed in de- 
tail. Traditional density functional theory (DFT) calcu- 
lations, based on local or semi- local approximations for 
the functional, provide detailed infor mation a bout the 
variations in atomic-adsorbate energies 22 23 24 but can- 
not account for dispersive interactions! 21 * 25 ^ Thus, tra- 
ditional DFT calculations have only a limited ability to 
compare binding energies, diffusion and rotation barri- 
ers for the molecular adsorption and does not permit an 
analysis of the adsorbate self-organization. In particu- 
lar, the generalized gradient approximation (GGA) has 
proved inadequate for Bz on Cu(lll). 27 In this paper 
we report vdW-DF calculations of the subtle variation in 
the adsorption energy of individual molecules with bond- 
ing site, of the intermolecular vdW binding, and of the 
energy barriers for adsorbate dynamics at different ad- 
sorbate condensations. 

Furthermore, we use the vdW-DF results to analyze 
and possibly interpret the self-organization mechanisms 
observed for Bz adsorption on Cu( 111) .^Assuming that 
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FIG. 1: [Color online] Schematics of the eight high-symmetry 
configurations of adsorbed benzene (Bz) on Cu(lll). High- 
symmetry configurations are found when the center of the 
molecules align with the traditional characteristic atomic- 
adsorption sites of the Cu(lll) surface (four main panels). 
In the illustrated I ("in-phase") orientation, one of the two 
equilateral triangles formed by alternate C's of the ring sits 
above three Cu atoms in the surface for HCP and FCC; for 
BR and TOP, these 3 carbons can be similarly aligned via a 
translation. Orientation R30° differs from orientation I by a 
30° rotation of the Bz adsorbate, as shown in the top pair 
of sketches. The arrow in upper left main panel illustrates a 
possible diffusion path that connects the high-symmetry sites, 
running from an HCP-I site, over a bridge site BR-I, then over 
an FCC-I site, and finally to a TOP-I site. A similar path 
could occur for the R30° orientation. 



the admolecules form a close-packed (hexagonal) struc- 
ture, we determine its lattice constant from the mini- 
mum of the interaction potential. For the more compact 
C2 phase, our vdW-DF calculations of the direct vdW 
binding indicate a preference for a lattice constant of 
6.9 A. Given the tendency of the vdW -DF m ethod to 
slightly overestimate binding separations, 1 21 1 28 1 this result 
suggests that the C2 phase can be interpreted as a direct, 
vdW-bound molecular-lattice structure. Explicit vdW- 
DF predictions of the preferred Bz orientations and of 
rotational barriers may offer experimental tests of this in- 
terpretation. For the less dense CI phase we find that the 
lattice spacing lies beyond the range of substantial bond- 
ing via the direct intermolecular vdW binding. How- 



ever, use a Harris-functional scheme^U formally extends 
our DFT analysis and shows that the CI lattice spac- 
ing agrees excellently with the fir st en ergy minimum of 
the indirect electronic interactio n 6 * 7 * 8 ^ mediated by the 
Cu(lll) metallic surface-state (Shockley) band. Since 
we have explicitly verified that the Cu(lll) surface state 
survives the addition of Bz molecules, we find it plausi- 
ble that this indirect electronic interaction underpins the 
self-organization in the CI phase. 

This paper has the following plan: in Sec. II we sum- 
marize the vdW-DF theory and formal Harris-functional 
evaluation that we use, along with various computational 
details. In Sec. Ill we present our results for adsorption 
energies, diffusion barriers, and barriers to rotation. In 
the next section we discuss these results, emphasizing 
possibilities for crystal organization. The final section 
contains our conclusions. 



II. DENSITY FUNCTIONAL THEORY 

A recent vdW density functional, vdW-DF) 19 | 2Q | has 
proved useful in describing a variety of sparse matter pi 
i.e. systems which contain important regions of low 
electron density where dispersive interactions contribute 
significantly to the materials' coherence. The vdW-DF 
method provides a transferable account of general materi- 
als and bindings, ranging, e.g., from the all-covalent C-C 
bonding within a graphene sheet to the interlayer dis- 
persive bonding in graphite and graphite intercalates.^ 
DFT calculations using the vdW-DF has successfully 
characterized the adsorption of Bz on the semi-metal 
graphene^ and on the semiconductors silicon 32 and 
MoS 2 . 21 A vdW-DF study also documents that vdW in- 
teractions (nonlocal correlations) significantly affect the 
anchoring of phenol on an oxide There are fewer ap- 
plications for the interaction between organic molecules 
and metal surfaces, but very recent examples include the 
adsorption on Cu(110) of triophene by Sony et alM^and 
of Bz and other organic molecules by Atodiresei et alP^ 

We describe the physisorption and two-dimensional 
assembly of Bz on Cu(lll) within DFT by combining 
large-scale vdW-DF calculations and a formal Harris- 
functional evaluation^ for: (a) adsorption of isolated 
Bz molecules and (b) direct and indirect (substrate- 
mediated) interactions between Bz adsorbates. Brute- 
force vdW-DF calculations treating the physisorption 
and the inter-adsorbate interactions in one giant unit 
cell are not computationally tractable for Bz adsorbates 
on Cu(lll). Instead, we apply large-scale vdW-DF cal- 
culations to characterize the molecular adsorption both 
in the case of a frozen surface and in the presence of 
adsorption-induced surface relaxations. A separate set 
of vdW-DF calculations investigate the direct molecu- 
lar interaction between pairs and in a regular three- 
dimensional lattice by ignoring the coupling to the sub- 
strate and by focusing strictly on the two-dimensional 
arrangement of adsorbates. Finally, we use a general, an- 
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alytical Harris- functional evaluation^!! for a nonpertur- 
bative treatment of the indirect electronic interaction be- 
tween adsorbates to formally complete an approximative 
DFT characterization of the adsorbate interaction and of 
the two-dimensional self-organization of Bz on Cu(lll). 
This indirect electronic interaction is oscillatory and very 
long range; it arises from interference of Friedel oscilla- 
tions produced by the strong adsorbate-induced scatter- 
ing of surface-state electrons! 5 1 8 1 36 1 For atomic adsorbates, 
the indirect electronic interaction can be described in tra- 
ditional DFT implementations using very large scale unit 
cells; 10 the present vdW-DF study must, in practice, rely 
on an implicit parameterization given by the detailed na- 
ture of the adsorbate-induced scattering! 7 * 8 * 37 ! 

Our splitting of the physisorption and assembly prob- 
lems permits us to provide explicit vdW-DF determi- 
nations for some properties but can, naturally, only 
serve as an approximate characterization of other prop- 
erties. Our description of the combined interactions of 
the adsorbates ignores, for example, the contributions to 
the direct vdW binding between molecules from small, 
adsorption-induced changes in the electron density distri- 
bution. Such effects should be present even if our use of a 
Bader analysis^* 3 - 9 -! documents that the total adsorbate- 
to-substrate charge transfer vanishes. Also, the precise 
nature of the surface-state scattering induced by Bz ad- 
sorption is not known. Lacking a detailed experimental 
characterization, we can only apply our formal Harris- 
functional evaluation to provide a qualitative d escri ption 
of the long-range indirect electronic interaction 5 * 6 * 7 * 8 * be- 
tween adsorbates. Nevertheless, we argue that the main 
qualitative results, possible interpretations of the CI and 
C2 phases, are robust. Moreover, the separation of the 
general physisorption and assembly problems should not 
significantly affect the results for the dilute-phase adsorp- 
tion energies, the minute surface corrugation, and the en- 
hancement of the C2 formation energy by direct Bz-Bz 
vdW bonding. 



A. van der Waals density functional calculations 

Since details of the vdW-DF method and of the func- 
tional design are given in Refs. 19 20 21 28, we here 
only summarize the calculational approach. For every 
atomic configuration in a unit cell Vb, we obtain the self- 
consistent charge with traditional DFT, typically as pa- 
rameterized in PBE-GGA. 40 From the known electron- 
density variation, we also calculate a local wavevector 
qo(r) which characterizes the energy density in a GGApHl 
Next, in a post-processing phase, we calculate a nonlocal- 
correlation terrrP^ 



Ef[n} = \ [ dr [ dr'n(v)<j>(r,r')n(v' 



) (i) 



as a spatial double integration over an interaction ker- 
nel 0(r, r'), where the latter integration is only over 



the unit cell Vq. The kernel can be conveniently eval- 
uated from a tabulated, universal form by introducing 
the two dimensionless variables S = (go + #o)l r — r/ | an( ^ 
£ = (go — Qo)/(qo + Qo)- The nonlocal correlation term 
E^ x [n) describes the interactions mediated via coupling 
to the electrodynamical field and allows vdW-DF to de- 
scribe sparse-matter behavior. We also evaluate a (semi- 
local) part E$ which is obtained from the PBE-GGA en- 
ergy by adjusting the exchange to the revPBE 41 form 
and the correlation to the value in in the local-density 
approximation!^ This term retains information about, 
for example, the kinetic-energy repulsion that adds to 
the surface corrugation? 3 -* Finally, we determine the to- 
tal vdW-DF energy as the sum 
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in a scheme that is transferable and capable of de- 
scribing inter action s for both high and low electron 
concentrations!^ 1 ^* 

In many respect the present vdW-DF calculations are 
similar to those presented in Refs. 1301311321331431441 The 
physisorption system requires us to resolve very small 
energy differences (e.g., the variation in adsorption en- 
ergy with site) and to ensure accurate and well- converged 
evaluations of changes in the nonlocal correlation E™ 1 [n] . 
Since Eq. involves a real-space integral, the evaluation 
of E^[n] depends slightly on the grid-samplin g beca use 
of the rapid variation in the kernel at small S ) 19 * 43 * We 
employ a number of additional calculational steps to en- 
hance accuracy and test the convergence in the evaluation 
ofEq. 0: 

We characterize a general adsorbate system L A by the 
Bz configuration a and the associated substrate config- 
uration S a (which may be the frozen free surface or the 
surface specified by adsorbate-induced relaxations as ob- 
tained by calculations and implementations of GGA-PBE 
forces). For a given atomic configuration A, the nonlo- 
cal correlation energy will be dominated by contri- 
butions from core regions at or close to the atomic po- 
sitions; this causes a very small relative error 5E^ l (A) 
that varies when atoms are displaced relative to the den- 
sity grid. To compute adsorption energies and compare 
two systems, A and B, the relevant nonlocal-correlation 
energy difference is instead specified by the smaller E^ 1 
contributions that arise outside the atomic cores. How- 
ever, since the two systems have different atomic config- 
urations (and hence different alignment with the density 
grid), the small relative errors [5Ef(A) and 5E^ l (B)} 
need not always cancel and could affect the the accuracy 
for systems with minute energy differences. In this study 
we therefore avoid direct evaluations of the nonlocal- 
correlation energy differences E^ l (A) — E^ l (B) and ap- 
ply, instead, the procedure described in Refs. 33 431 In 
practice, we determine E^ l (A) — E^(B) in a sequence of 
steps where we separately calculate the E^ 1 changes as- 
sociated with removing the Bz in configuration a (b) to 
infinity, always keeping the atomic positions frozen and 
fixed relative to the density grid. In the vdW-DF study 
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of the effects of adsorbate-induced relaxations of the sub- 
strate, we follow these steps by a GGA-PBE calculation 
of the energy difference between S a and S5. A previ- 
ous study^ documents that a grid spacing of 0.15 A is 
sufficient to achieve convergence of the vdW-DF binding 
energy in adsorption problems when combined with the 
above-mentioned scheme. 

The calculations of E^ l [n] via Eq. |l]) are accelerated 
by introducing a cut-off radius R cut while carefully test- 
ing the convergence of the spatial double integration. 
Contributions of the long-range dispersive interactions 
to E^ l [n] fall off slower with distance |r — r'| for the 
adsorbate-slab interactions than for the Bz-Bz interac- 
tions. The r 7 integration in Eq. ([I]) exceeds the unit cell 
but the integrated weight of contributions outside Rcut 
still decreases fairly rapidly with the value of R C ut^^ 
We have explicitly tested that the choice R cut = 23 A 
ensures that E£ l converges at the sub-meV level. 



B. Long-range indirect electronic interactions 
between adsorbates: nonperturbative s-wave 
treatment in the asymptotic regime 

To complete our DFT analysis we must explicitly add 
consistent estimates for the indirect, substrate- mediated 
interactions. Benzene adsorbates not only interact di- 
rectly through vdW forces, but also through forces me- 
diated by the Cu(lll) substrate. Elastic-deformation 
effects,^ dipo lar couplings,^ and indire ct e lectronic 
interactions^^ all exhibit a power-law deca)^^ and can 
affect the adsorbate binding at distances d > 8 A where 
(we find) the direct adsorbate-adsorbate vdW interaction 
is no longer substantial. Since we calculate the direct Bz- 
Bz vdW interaction in a separate vdW-DF study that 
ignores the substrate, we must separately analyze the set 
of indirect interactions and add the relevant DFT-based 
terms. 

Of particular interest is the surface-state-mediated in- 
direct interaction. This metallic surface state is respon- 
sible for the dramatic images obtained with scanning- 
tunneling microscopy (STM)™ Fig. H shows that there 
is notable coupling between the Bz adsorbate and the 
surface-state wavefunction even in a dilute-coverage 
regime. This conclusion is expected from the analogous 
case of atomic adsorption^ on Cu(lll) 5 where the mea- 
sured adsorption-induced phase shifts are known to be 
large - 9 = 9f~tt/2 — indicating that a perturba- 
tive treatment of the surface-state-mediated adsorbate- 
pair interaction is inadequate. Our vdW-DF calculations 
of Bz physisorption suggest that the dipolar and elastic 
interactions between Bz adsorbates should be insignifi- 
cant: a Bader analysis^^ documents vanishing electron 
transfer and we find only very small adsorbate-induced 
relaxations, ~ 30 meV, which do not change the relative 
stability of the various sites. However, Fig. [2] shows that 
the metallic surface state survives the Bz adsorption and 
therefore can mediate an interaction between adsorbates. 
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FIG. 2: Robustness of Cu(lll) surface- state formation to Bz 
adsorption and evidence of adsorbate-induced surface-state 
scattering presented in the dilute limit with one Bz molecule 
on a frozen 3x3 Cu(lll) supercell. The figure shows the vari- 
ation of the in-plane average of the surface-state wave func- 
tions squared, |^| 2 = A -1 J 2Dcell dA iVWf-st. | 2 (evaluated at 
the f point) along the surface-normal coordinate z for our 
choice of a 6-layer copper slab (with top and bottom Cu layers 
identified by light dashed lines) . The physisorption position of 
the Bz adsorbates is indicated by the black dashed line (near 
z — 19A). Many layers are necessary to fully describe the 
surface- state properties and, in a slab geometry, to decouple 
the surface state on the top and on bottom sides. Our choice 
of 6 copper layers proves to be the smallest that exhibit clear 
signs of a Shockley surface state — an exponential decay of 
the surface state into the bulk and a large spill-out far be- 
yond the outer layers. The approximative symmetry of the 
top-side (solid curve) and bottom-side (dashed curve) surface 
states implies that the surface states survive the physisorption 
in the dilute limit. At the same time, the oscillatory behavior 
evident at the position of the adsorbate shows that the Bz 
interacts with surface-state electrons. 

Fig. [2] also demonstrates that Bz affects the surface-state 
scattering and produ ces, in turn, a long-range indirect 
electronic interaction! 5 * 6 * 7 * 8 ^ The minor adsorbate-induced 
relaxations will only enhance the surface-state scattering 
by increasing the coupling between the surface-state and 
evanescent bulk-like states around the adsorption site. 

Unlike the other adsorbate-adsorbate interactions, the 
indirect electronic interaction is oscillatory. At short 
range, it depends on all the occupied states but asymp- 
totically the behavior is dominated by the q-state at the 
Fermi level whose velocity is parallel to the separation 
vector d between the adsorbates P While this interaction 
is typically anisotropic, the metallic (Shockley) surface- 
state band on Cu(lll) is circularly symmetric, so the in- 
teraction becomes isotropic, and the relevant is paral- 
lel to d. These conditions lead to the perturbative treat- 
ment well known as the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction* 48 * 49 ! which has recently been ap- 
plied to actively- investigated low-dimensional systems!^ 
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A significant effect of mediation by metallic surface states 
is that the envelope decays exceptionally slowly, ex d~ 2 
rather than the u sual oc d~ 5 for mediation by bulk states 
at metal surfaces! 5 * 6 * 7 * 8 * 

The primary modification due to a nonperturbative 
treatment is the addition of a phase shift in the argu- 
ment of the sinusoidal function in the RKKY expression 
(and for Friedel oscillations in general). To explicate this 
result, we extend our regular vdW-DF calculations with 
a Harris- functional scheme, 29 which formally constitutes 
an asymptotic DFT evaluation^ of the indirect electronic 
interactions between Bz adsorbates. This choice permits 
a consistent description of the nonperturbative scattering 
and interaction effects arising with general (complex) s- 
wave surface-state phase shifts. 7 The complex phase shift 
reflects transitions both among surface-state wavevec- 
tors and into evanescent bulk states. A simple s-wave 
characterization, 8 with simultaneous descriptions of both 
scattering channels, has proven reasonably successful in 
estimating the resulting indirect-electronic interaction for 
a pair of atomic adsorbates^ the s-wave characterization 
should also be applicable for Bz on Cu(lll) because the 
coupling of Bz to the substrate occurs in a relatively con- 
fined area. 

The s-wave formulation defines an effective Fermi-level 
T-matrix to(ep) oc [exp(i2c)o) — 1] given by a complex 
phase-shift So = 5 + i5" and formally expressing the 
probability for transition entirely within the surface- 
state channel. Unlike for a strictly 2D scattering sys- 
tem, however, the intra-surface-state (ISS) scattering is 
given both an ISS phase shift 5 and an ISS reflectance 
r = exp(— 25"). The latter corresponds to a bulk-state 
absorption probability 1 — r, which cannot vanish since 
the adsorbates cause surface-state electrons to scatter 
into bulk states. The nonperturbative analysis leads to 
an effective Fermi-level phase shift 6 which can be ob- 
tained frorrP 
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where the expansion is around either endpoint or 7r/2, 
with tp = S or (7r/2) — 5, with the upper or lower sign 
applying, respectively. For small r, G(r, 5) is nearly sym- 
metric about S = 7r/4, descending linearly (from value 
7r/2), with a shallow slope, from either endpoint. As r 
increases, the magnitudes of the linear slopes increase 
asymmetrically, with a shift of the minimum 5 m i n to pro- 
gressively smaller values, while the minimum value of O, 
G(r, £min) 5 decreases continuously from tt/2 to 0; explic- 
itly, 



x 1 -i 

Omin = - COS T, 



@(Mmin) = tan 



(4) 



The upper panel of Fig. [3] depicts the relation between G 
and 6 for three representative r values. 
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FIG. 3: [Color online] Dependence of effective interaction 
phase shift G (top panel) and scaled interaction strength 
A of indirect electronic adsorbate-pair interaction [Eq. |5])] 
with detailed nature of adsorbate-induced s-wave scattering 
of surface- state electrons. The variation is expressed by the 
phase-shift for intra-surface-state transitions S and the intra- 
surface-state reflectance r (where 1 — r describe the probabil- 
ity that the adsorbate causes a scattering from a surface state 
into an evanescent bulk- like state). Analogy with the case of 
interactions between Cu(lll) adatoms suggests that it is rea- 
sonable to expect both 1 — r > and 5 > 0. While we do 
not have explicit determinations of the relevant parameters 
for Bz on Cu(lll), the panels suggest that values, G > tt/3 
and A > 0.1, can be considered typical for indirect electronic 
adsorbate pair interactions [valid for the present case of Bz 
on Cu(lll) as well as for Cu adatoms on Cu(lll/^. 



We determine E*™? r (d), where the superscript de- 
notes that this is the asymptotic form of the indirect- 
electronic adsorbate-pair interactions, computed with the 
nonperturbative formal Harris-functional for s- wave scat- 
tering, through the analytical evaluation^ 



/4e F \ sm(2q F d + 2S(r,S)) 



{q F df 



(5) 



Here 



A(r, S) = [(1 - r) 2 /4 + r s'm 2 (2S)} (6) 
represents a scaled, dimensionless interaction strength, 
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< A < 1. For Cu(lll) the wavevector qp of the surface 

state is 0.22 A , and the Fermi energy ep is 0.39 eVP 
The explicit values of S and of r, or equivalent ly 9 and 
A, depend on the details of the Bz-induced scattering, 
for which we have no experimental characterization; we 
stress that, in the present s-wave analysis, the effective in- 
teraction phase is the same O that characterizes the STM 
measurements of the local density-of-state variationP^ 

Fig. [3] displays the possible variation in interaction 
phase 6 (upper panel) and scaled interaction strength 
A (lower panel) with the s-wave modeling parameters 5 
and r. The perturbative regime of RKKY corresponds to 
the limit (r —> 1, 5 —> 0), while r — > corresponds to the 
case of completely absorbing scattering (all surface-state 
electrons scatter into the bulk), for which Eqs. Q and (|6| 
indicate O = tt/2 and A = 0.25, respectively. In making 
these plots we have implicitly assumed a positive inter- 
action phase shift 6 > (as was measured for atomic 
adsorbates^). We expect Bz adsorbates to produce sig- 
nificant coupling between surface and bulk states, r / 1 
and finite ISS scattering 5^0. The lower panel of Fig. [3] 
shows a contour plot of the scaled interaction strength A; 
even moderate values of the ISS scattering phase-shift S 
and of the bulk-state absorption, 1 — r, produces signifi- 
cant values of the interacting strength. Generic values for 
the ISS scattering parameters (5 and r) produce an in- 
teraction strength and phase which are comparable to or 
larger than those characterizing interactions between Cu 
adatoms on Cu(lll)P A Cu /Cu = 0-08 an d ©Cu/Cu = ?r/3. 

It is interesting to compare the indirect electronic in- 
teractions between pairs of Bz and pairs of Cu adatoms 
on Cu(lll). In both cases is there is essentially no charge 
transfer. The Friedel sum rule 51 then tells us^ that the 
coupling between surface and bulk-like states must be 
the primary cause of the Friedel oscillations. On the 
one hand, since Bz is physisorbed to Cu(lll), its im- 
pact should be smaller than chemisorbed Cu adatoms, 
with a smaller adsorbate-induced coupling between the 
surface and bulk statesP^ On the other hand, Bz has a 
quadrupole moment, defined by a total of 6 C-H bonds 
with a nonzero dipolar moment that will scatter the 
surface-state electrons. In a simple tight-binding model 
of weak adsorption, the indirect interaction varies as the 
square of the adsorption energy.^ It is plausible that scat- 
tering off each carbon atom in Bz produces an indirect 
electronic pair interaction with a nonvanishing interact- 
ing strength and phase, so that there is a significant 
(though likely not 36-fold) enhancement of the individ- 
ual interactions. One might instead view the adsorption 
as due to a molecular orbital of each cyclic group, im- 
plying that the indirect interactions of anthracene pairs 
might be 3 - 9 times those of benzene at the same separa- 
tion. Lukas et alP^did find the intriguing result that the 
binding energies of C4 n +2H2 n +4 on Cu(lll) for n=0,l,2,3 
(acetylene, benzene, naphthalene, anthracene) are ap- 
proximately collinear as a function of n, with the binding 
energy of anthracene about 5 times that of acetylene and 
somewhat more than twice that of benzene. In summary, 



it is not clear how several adsorbates within the screening 
length of the bulk electrons couple to the surface state 
compared to just a single adatom. This issue begs further 
investigation. 



C. Adsorption and surface assembly energies 

The adsorption energy is given by the difference be- 
tween the total energy for a given configuration and the 
total energy of the Bz and the surface in isolation 



K 



vdWDF 

ads 



E 



vdWDF 



(h) - E 



vdWDF 



(h - 00). (7) 



(In practice, due to the finite size of the unit cell, each 
term on the RHS of Eq. (7) contains ~30 meV per 
Bz due to lateral vdW interactions, which essentially 
cancel after the subtraction.) The full formation en- 
ergy when including the energy gained by forming an 
overlayer is given by the sum of the adsorption energy 
and the two-dimensional (2D) assembly energy Ef orm = 
E a ds + -E^Dassembiy • I n this paper, we approximate the 
assembly energy of a 2D overlayer by the pairwise sum 
of Bz-Bz interactions: 



T^app _ 
■ C/ 2Dassembly ~ 2 / > ^pair I* tj 



( r »i) • 



(8) 



This corresponds to neglecting higher-order (multi- 
molecule) Bz-Bz interactions both for the van-der Waals 
and indirect electronic interactions. A test calculation 
with a moderately dense benzene overlayer indicates that 
changes in the vdW-DF result can alter the formation en- 
ergy by a few tens of meV. For indirect electronic inter- 
actions at large separations, the contribution of trio (3- 
adatom, non-pairwise) interactions decays only slightly 
faster than pair interactions and should be, very roughly, 
1/4 the magnitudeP^ While such higher-order interac- 
tions are needed to achieve high accuracy, they are not 
essential for our analysis. We further argue that to a good 
approximation elastic effects and surface corrugation can 
be neglected in computing inform • We thus calculate the 
direct component of the Bz-Bz interaction in isolation 
from the substrate: E^^Zoo- A 

gain this energy is 
the difference between nearby molecules and widely sep- 
arated ones: E;™ DF (d) = £ vdWDF (d) - £ vdWDF (d - 
oo). The interactions between molecules mediated by the 
surface, formally included in Eo, is reintroduced by in- 
cluding the leading-order asymptotic part of the surface- 
state-mediated indirect interaction as discussed in the 
previous section: 



pair 



? vdWDF 
^pair,fi— >oo 



T^iaris 
-^pair * 



(9) 



The asymptotic expression E^f Y is expected to provide a 
good approximation to the full pair indirect interaction 
in the regime where the latter dominates the direct vdW 
interaction. 
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D. Choice of unit cell, computational details 

The self-consistent charge density is found using the 
plane- wave and ultrasoft-pseudopotentials code DACAPC 52 
with GGA in the PBE flavor. For the physisorption study 
on a 3x3 slab unit cell with 6 copper layers, we use a 
Monkhorst-Pack (3,3,1) k-point sampling. The study of 
direct vdW-interaction between Bz pairs in the dilute 
case was done in a very large unit retaining only the T 
point. The plane- wave energy cutoff was set higher than 
500eV to ensure that the spatial gridding was consistently 
better than 0.14 A. 

We place a Bz molecule on a 3x3 Cu(lll) supercell 
with a slab containing about 6 copper layers. The choice 
of in-surface dimensions corresponds to a Bz-to-Bz sep- 
aration of a ce ii = 7.74 A. Our system is large enough to 
extract the essential physics but not so large as to pre- 
clude doing so in a timely fashion. Our vdW-DF calcu- 
lations are based on a periodic unit cell and the choice of 
supercell does leave a small (but fixed) residual molecule 
coupling arising from both direct inter-molecular vdW 
interactions and from substrate-mediated interactions. 
The vdW-DF study of the direct vdW-interaction be- 
tween Bz pairs indicates that the effect of direct vdW 
bonding is small (3 times ~ 10 meV); this energy is as- 
sumed to cancel in the subtraction in Eq. (7). However, 
the long-ranged (and oscillatory) nature of the indirect 
electronic interaction motivates a separate discussion as 
do the convergence of physisorption energies with slab 
thickness. 

We note that the effects of indirect electronic 
adsorbate-adsorbate interaction will be contained in 
very-large-scale DFT calculations^ and that vdW-DF 
would formally contain such a contribution, primarily in 
the semi- local term Eq. However, our specific choice of 
supercell with 6 copper layers does not allow more than a 
qualitative discussion of the surface state, Fig. 2, and it is 
not possible to use our vdW-DF calculations for a quan- 
titative analysis. Instead we note, based on an assump- 
tion of completely absorbing surface-state scattering in 
the asymptotic limit, that the magnitude of the fixed 
contribution from indirect electronic interactions (arising 
due to our use of an implicit periodicity in our vdW-DF 
unit cell) can be approximated as 3x2=6 meV (accord- 
ing to Eq. (5), which is not fully valid at this not-quite- 
asymptotic separation). This value is much smaller than 
the calculated vdW-DF value for the Bz physisorption en- 
ergies and we choose not to correct for the fixed effect of 
implicit indirect electronic interaction in presenting our 
vdW-DF results for the isolated- molecule adsorption. It 
is safe to assume that these inter-unit-cell contributions 
will not modify the vdW-DF results of extremely low dif- 
fusion barriers and essentially corrugation-free dynamics 
of individual Bz adsorbates. 

A partially related issue is convergence of vdW-DF en- 
ergy differences with slab thickness. We find the variation 
of the physisorption energy with adsorption site to be 
very small, a finding that motivates us to enhance the ac- 



curacy in the vdW-DF calculations. To test convergence 
with respect to slab thickness, we compared the results of 
vdW-DF calculations for two different adsorbate geome- 
tries. Specifically, we have determined how many Cu lay- 
ers n are needed to accurately and reproducibly predict 
energy differences for adsorption configurations FCC-I 

QT1 j rpn pono. A pvdWDF _ pvdWDF rvdWDF 

and r^O-rUU . l\^ n - E FCC _ ln - ^ FC C-R30,n- 

Here, Ai^ dWDF changed from 10% when comparing 4 
and 5 layers to 2% when comparing 6 and 7 layers. We 
thus conclude that 6 copper layers give satisfactory con- 
vergence. The somewhat slow convergence might be re- 
lated to the need for several Cu layers before a clear sig- 
nature of the Shockley surface begins to emerge. 

III. SINGLE MOLECULAR ADSORPTION AND 
DIFFUSION 

A. Physisorption energy 

Fig. [4] displays the potential energy landscape, found 
by performing a potential energy search along an axis 
of high-symmetry sites, as shown in Fig. [TJ for both BZ 
orientations (configurations) I and R30°, as defined in 
the caption of that figure. We allow the molecule to 
adjust its height at each step, but do not allow further 
substrate relaxation. The spline- interpolated result for 
the I configurations are shown in the upper panel, while 
the lower panel shows the binding curve for four of the 
high-symmetry sites. All I configurations have a lower 
energy than all R30° configurations, and the line between 
HCP-I, BR-I and FCC-I sites is found to be almost de- 
generate, with a predicted adsorption energy of 0.51 eV, 
where FCC-I is slightly more favorable than the two other 
sites 

The good agreement with the experimental binding en- 
ergy 0.6 e V 12 | 13 [ further demonstrates that the use of the 
vdW-DF offers a viable way to describe the binding be- 
tween organic molecules and noble metal surfaces, even 
in the presence of a metallic surface state. The accuracy 
of using vdW-DF for this system is similar to its per- 
formance in general. The accuracy is also close to that 
found for n-butane on Cu(lll)P^ 

Our results clearly demonstrate the non-local nature of 
the binding. The binding separation of 3.75 A is a typical 
separation distance for physisorption; furthermore, with 
a Bader analysis^ 3 -^ we found vanishing charge transfer 
between the surface and the Bz molecule. 

In the results presented so far, we use a frozen sub- 
strate, thereby neglecting elastic effects. To assess the 
importance of these effects, we perform for every high- 
symmetry site a set of separate calculations, in which we 
let the top 4 layers of the surface relax and fix the coordi- 
nates of the benzene atoms in positions ranging from 0.1 
above to 0.1 A below the optimal separations to the un- 
relaxed surfaces. We use PBE-GGA for this relaxation, 
while in the post-processing phase, the relaxation energy 
is estimated from vdW-DF. 
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FIG. 4: [Color online] Potential-energy curves for Bz on 
Cu(lll), for computing diffusion barrier. For the contour plot 
shown in the upper panel, the molecule is allowed to change 
its separation from the surface at each position. The vertical 
axis gives the separation from the surface and the horizon- 
tal axis shows the translation along the arrow indicated in 
Fig. ^ The bar to the right shows the energy scale in meV 
for the equipotential lines. This plot was smoothed using 
spline-interpolation. The lower panel shows the potential en- 
ergy curves for Bz adsorption in the four high-symmetry sites 
with I orientation. The upper lower (upper) solid curve is for 
the FCC-I (TOP-I) site, the dashed (dotted) is for the HCP-I 
(BR-I) site. The insert in the lower panel shows the nonlocal 
part of the potential energy curves, which are very similar, 
illustrating that the energy difference between the sites stems 
primarily from .Eo[n]. 



Table [T] give the adsorption energy of the eight high- 
symmetry sites, for a frozen substrate, a substrate al- 
lowed to relaxed and for the purely PW91-GGA calcu- 
lation reported in Ref. [27j It shows that allowing re- 
laxation lowers the adsorption energy within vdW-DF 
by between 29 and 34 meV, getting even closer to the 
experimental value of 0.6 eV. Allowing relaxation gives 
only minor corrections to the relative favorability of the 



TABLE I: Interaction energies for adsorption of benzene 
(CeHe) at various positions on Cu(lll) in order of decreasing 
stability, with vdW-DF but frozen substrate, with vdW-DF 
and relaxation, and without vdW-DF [from Ref. 127]. 



Site 


^vdW-D F[meV] 


7-wclW-DF r T7i 
^relaxads. [ meV l 


^ds A [meV] 


FCC-I 


-515 


-548 


-23 


HCP-I 


-514 


-548 


-26 


BR-I 


-514 


-547 


-16 


TOP-I 


-503 


-533 


44 


HCP-R30 


-499 


-533 


-18 


FCC-R30 


-498 


-527 


-17 


BR-R30 


-497 


-529 


-21 


TOP-R30 


-490 


-520 
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different adsorption sites. Note, though, that the calcu- 
lations with relaxation are more sensitive to the shifting 
from the PBE to the revPBE flavor of exchange in the 
post-processing phase, and so are converged only at the 
2-3 meV level rather than the sub-meV convergence of 
the calculations without adsorption-induced relaxation. 
Since the better convergence gives us greater confidence 
in the computations with a frozen substrate, we base our 
investigation of the surface-corrugation results on them. 

Finally, Table [I] shows dramatically the need for includ- 
ing vdW interactions to accurately describe this adsorp- 
tion system. Calculations which do not include vdW [in 
this case^, but use PW91-GGA done with the Vienna 
Ab-initio Simulation Package (VASP)] underestimate the 
adsorption energy by a factor of around 20. For the two 
TOP configurations, the benzene is not even bound. The 
near degeneracy of the three lowest configurations is lost, 
and the general ordering is considerably different. The 
theme of the overall favoring of I-configurations is absent. 
Thus, such calculations are inadequate for semiquantita- 
tive and perhaps even qualitative purposes. 



B. Surface corrugation, barrier-less diffusion on a 
honeycomb lattice 

Fig. [4] presents our vdW-DF results for diffusion bar- 
riers and shows that Bz adsorbates are expected to opti- 
mize their positions relative to their mutual interactions, 
with the substrate corrugation potential playing a rel- 
atively minor role. Table [I] presents a overview of the 
surface corrugation as calculated in vdW-DF and GGA. 
We note that the surface corrugation (and hence diffu- 
sion barriers) is at least an order of magnitude larger 
in GGA but we do not believe that the GGA finding, 
top sites being unstable towards desorption, is physical. 
The vdW-DF calculations predict a much smaller cor- 
rugation: a low-temperature preference for placing Bz 
adsorbates around the nearly-degenerate FCC-I, HCP-I, 
and BR-I sites with a relative cost of around 11 meV for 
placing Bz on a TOP-I site. Our vdW-DF description 
may slightly underestimate the corrugation because it is 
likely to weakly overestimate the physisorption height h. 
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The TOP-I versus FCC-I/HCP-I/BR-I energy difference 
is in any case sufficient to affect the expected equilibrium 
distribution at 77 K. Our vdW-DF results shows that 
the Bz adsorbates at low temperatures will concentrate 
in a honeycomb pattern (physisorption valley) formed by 
linking the FCC-I, HCP-I, and BR-I sites. 

At the same time, we find that the Bz adsorbates ef- 
fectively experience no diffusion barriers and are free to 
respond to their mutual interaction. All of the diffu- 
sion barriers are tiny compared to those for chemisorbed 
atoms; for a Cu(lll) adatom the experimental and the- 
ory results for the diffusion barrier is 40 meVP^ 23 | 24 [ 
Even for a hypothetical diffusion path, HCP-I — > TOP- 
I — > FCC-I, the peak of the barrier (11 meV) is small 
enough to allow a diffusion rate larger than 1 s _1 at 
T > 4 K. The diffusion barrier for motion within the 
honeycomb pattern, along the path HCP-I — » BR-I — > 
FCC-I, is less than 1 meV and we stress that the Bz 
adsorbates can move continuously across the surface, us- 
ing this bend-like type of diffusion. Again, our vdW- 
DF calculations may slightly underestimate the diffusion 
barriers but it is clear that the Bz adsorbates are free 
to respond to their mutual interactions on the Cu(lll) 
surface. 

Despite such low barriers, there are many reports of 
lattice-gas phases of physisorbed atoms at submonolayer 
coverages. 54 For example, while there is an incommen- 
surate phase of Kr on graphite, there is also a widely- 
studied commensurate phase (^/3 x ^3)^30° with one in 
three sites occupied on the hexagonal grid formed by the 
centers of the honeycombs.^ The theoretical analyses^ 
of the disordering phase transition involves a lattice-gas 
model appropriate to adsorption on a grid of sites. Like- 
wise, we find such behavior below for the two condensed 
phases of Bz on Cu(lll). 



C. Rotational barriers 



Having established that the benzene can move almost 
freely on the surface of Cu(lll), we investigate whether 
it can also rotate freely. The symmetry of the I and 
R30° configurations and the monotonic increase in the 
barrier profile from HCP-I/FCC-I to TOP-I allow us to 
assume that the peak of the rotation barrier is given sim- 
ply by the energy difference between orientation I and 
orientation R30°. 

Table [TT| compares the energy differences between Bz 
adsorbates in orientation R30°and I, which in turn deter- 
mine the rotational barriers. The vdW-DF results show 
that individual molecules can easily rotate, even though 
the associated activation energy is slightly higher than 
for diffusion. 



TABLE II: Binding difference between R30°-sites and I-sites, 
which also give the rotational barriers. The second column 
show the results for a frozen substrate, while the third column 
gives the results when relaxation is included. 



Site 


A£ a v d d s WDF [meV] 


A£ r v e 1^d F s[meV] 


HCP 


15 


15 


BR 


17 


18 


FCC 


17 


21 


TOP 


13 


13 



IV. SURFACE ASSEMBLY 

The vdW-DF findings of small diffusion barrier and 
lack of expected surface corrugation for individual Bz 
adsorbates support the picture of an essentially free 2D 
gas of Bz molecules in the dilute regi me. T his is what 
was found on Cu terraces at T = 77 K.™ Further, the 
finding that the location of Bz adsorbate to some extent 
decouple from the underlying copper lattice, leaves the 
stage open for other phenomena to dominate the inter- 
play between molecules at higher coverages and deter- 
mine the 2D surface assembly of the Bz adsorbates. 

Dougherty et alP^ observed the coexistence of a 
high density (2.54nm -2 ) phase C2 and a low density 
(l.lOnm -2 ) phase CI, both forming hexagonal patterns. 
We seek insight and provide an analysis in which we first 
calculate the direct and indirect Bz-Bz interactions and 
then combine these into an effective interaction poten- 
tial. We finally use this interaction potential to charac- 
terize the high-density phase and discuss the possibility 
of surface-state-mediated forces being responsible for the 
structure of the low density phase. 

Table |III| summarizes the properties of the three Bz- 
on-Cu(lll) regimes. These results can be compared with 
measurements to test the accuracy and predictive power 
of our vdW-DF calculations and the analysis of the indi- 
rect elect roni c interaction mediated by the Cu(lll) sur- 
face state.™ 



A. Results for pair interactions 

There are two important forces between the molecules 
for the interacting system formed when multiple Bz ad- 
sorbates assemble on Cu(lll). First is the direct vdW 
interaction between the Bz molecules; this interaction is 
medium range. Second is the indirect electronic Bz-Bz in- 
teraction, mediated by the surface-state electrons, which 
is long-range and oscillatory! 5 * 6 * 7 * In this subsection, we 
address the influence of these effects for pair interactions 
and construct an effective pair potential which can be 
used to discuss island formation and phase condensation 
of Bz on Cu(lll)P 

Fig. [5] shows the vdW-DF binding curve of Bz pairs, 
^pairT3oo(^) f° r different high-symmetry configurations 
identified in the set of inserts and calculated in vdW-DF 
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TABLE III: Energy and structure of the different phases. The 
formation energies (calculated from the expression for inform 
after Eq. (7)) are listed both as obtained for a frozen surface 
and (in parentheses) in the presence of adsorbate-induced sur- 
face relaxations. The rotational barriers are given by the site- 
orientation dependency of the benzene molecule. For CI and 
C2 data the molecules are assumed to be in the I orientation, 
with the effect of indirect electronic interactions calculated 
for A=0.25 and G = tv/2 (corresponding to a choice of com- 
pletely absorbing scattering). The square brackets give the 
values with A — 0. The rotational excitation is calculated 
under the assumption that the I- R30° interaction is limited to 
neighboring molecules. 



C2 



CI 



Dilute 



Formation energy [meV] 
Exp't adsorp'n energy [meV] 

Lattice parameter [A] 
Exp't lattice parameter [A] 

Rotation barrier [meV] 
Rotation excitation [meV] 
2D-assembly energy [meV] 



618 (-) 539 (572) 515 (548) 
600 a 



6.9 
6.8 b 

~ 9 
^79 [97] 



10.2 
10.2 b 

-17 
-24 



17 



Based on Ref. H b Based on Ref. fT5l 




FIG. 6: [Color online] Effective binding potential of a Bz pair 
from the combinations of direct vdW and indirect electronic 
(surface-state mediated) Bz-Bz interactions, assuming favor- 
able orientations of Bz molecules. The upper dashed curve 
shows the results for the values of A = 0.08 and G = 0.37T, 
while the lower solid curve shows the results for A — 0.25 
and G = tt/2. For d > 8.1 A (indicated by the dotted black 
line), the vdW interaction is negligible, and the displayed in- 
teraction is given by the formal Harris-functional evaluation, 
Eq. |5|. 
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FIG. 5: [Color online] Binding curve of two Bz molecules in 
different high-symmetry configurations identified by the set 
of inserts and calculated in vdW-DF in isolation from the 
surface. The two Bz's lies entirely within a common plane 
and the abscissa gives the distance between the center of two 
Bz's. We label the configurations shown in the lower left, 
lower right and upper right corner as acene-like (A), mixed(M) 
and phenyl-like(P). The most favorable A configuration (solid 
curve) has an optimal binding distance of 6.9 A, and corre- 
sponding vdW-DF binding energy is 32 meV. 



ignoring contact with the surface. The two Bz molecules 
lies entirely within a common plane. For lack of bet- 
ter notation in the litterature, we label them acene-like 
(A), phenyl-like (P) and mixed (M). In the A (P) con- 
figuration the displacement vector from one Bz to the 
other goes perpendicularly through an edge (a vertex) of 
both carbon hexagons. In the M configuration the dis- 



placement vector goes through a vertex of one Bz and an 
edge of the other. The binding energy is largest for the 
acene-like configuration, closely followed by the mixed, 
with the minimum of the bonding curve at —32 meV. In 
the P configurations bonding is much weaker and at a 
larger separation. The direct vdW binding clearly has a 
larger magnitude than the barrier profile for the adsorbed 
molecules. Thus, the Bz-Bz pair interactions should dom- 
inate over surface corrugation effects at smaller separa- 
tions (d smaller than ~ 8 A). Summing the contributions 
of pairwise vdW bindings in A configurations, we esti- 
mate the full vdW-assembly energy of an isolated sheet 
of Bz molecules in a hexagonal lattice to be 97 meV. 

On the Cu(lll) substrate, the I-configurations are op- 
timal for single-molecular adsorption and the A config- 
urations are optimal for pairs. It is thus relevant to 
discuss an effective adsorbate-pair interaction £? pa ir with 
implicit assumptions of favorable molecular orientations 
(with both Bz molecules in surface orientation I and with 
the direct vdW interaction specified by the A configura- 
tion). This effective (pair) potential determines adsor- 
bate organization into a well-ordered triangular lattice 
structures at very low temperatures. 

Fig. [6] shows our result for the effective (favorable- 
condition) potential ^ pa ir with a binding curve that 
includes the asymptotic surface-mediated interactions 
given by Eq. ([5| . Since we lack explicit information about 
the detailed nature of surface-state scattering by Bz ad- 
sorbate, we must make assumptions for it by analogy 
with the case of atomic adsorption on Cu(lll)P 

We give the results for the case of a completely ab- 
sorbing scatterer with A = 0.25 and = 7r/2, and, for 
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FIG. 7: [Color online] Schematic illustration of a rotational 
excitation in the hexagonal phase C2 and of the interlocking 
effects which (at low temperatures) will prevent neighboring 
adsorbates from rotating until the original rotational excita- 
tion relaxes. The figure shows two near-degenerate configu- 
rations of benzene on Cu(lll), in the left panel all molecules 
are in the ground-state packing with all adsorbates in the I- 
orientation, while in the right picture one of molecules has 
switched to the R30°-orientation. This switch corresponds 
to six adsorbate pairs changing to the M configuration. For 
pairs, this switch corresponds to six of them changing to the 
M configurations. While this rotational excitation does not 
cost much energy per se, it does block the rotation of all the 
neighboring molecules (marked in light blue [gray in print]) 
until the center molecules relaxes back to orientation I. 



comparison, the values fitted to the experimental scatter- 
ing of copper adsorbates 5 with A = 0.08 and = 0.37T. 
The black dotted vertical line in the figure separates the 
region in which van der Waals forces are important from 
the regions where the surface-mediated interactions dom- 
inate. For the latter region only the indirect-electronic 
interactions are included. For A = 0.25, the curve has 
two minima, the first stemming from direct intermolec- 
ular interactions, and the second, weaker one stemming 
from the sinusoidal form of the contribution For 
A — 0.08, the indirect-electronic interaction only creates 
a second, very weak minimum (1.5 meV) at 14 A, but 
do still create an energy barrier for the adsorbates try- 
ing to move from the CI phase and assemble into the 
vdW-bound C2 phase. 



B. Van der Waals lattice (C2) 

The optimal pair separation of 6.9 A, as shown in 
Figs. [5] and [6j illustrates that it can be favorable for is- 
lands to form at benzene separations somewhat shorter 
than that of the 3x3 unit cell: 3a ~ 7.74 A. The primary 
limitation for the density of these islands is the kinetic- 
energy repulsion between benzene molecules. Assuming 
that the molecules lie flat on the copper surface and form 
a hexagonal pattern with pairs of molecules in A config- 
urations, we estimate a density of 2.42 nm -2 , very close 
to the experimentally determined density. 

Fig. [7] illustrates a curious consequence for the adsor- 
bate dynamics in phase C2. We include it as the phe- 



nomenon is a consequence of the proposed interpretation 
of phase C2 as a planar hexagonal pattern of Bz adsor- 
bates, primarily organized in orientation I, with pairs of 
Bz's in the A configuration and assembled by direct vdW 
bonding; an experimental observation of this consequence 
would constitute a specific test of the proposed interpre- 
tation of phase C2 as a vdW-bound assembly. Based on 
vdW-DF calculations of the direct bonding, we find that 
the hexagonal pattern will permit a low-energy excitation 
with a rotation of a single Bz adsorbate from orientation 
I to R30°. This switch corresponds to six Bz pair configu- 
rations changing from A to M. However, as illustrated in 
the right panel, such a rotational excitation will (at low 
temperature) cause a blocking of rotations for all neigh- 
boring Bz adsorbates. This effect is due to the strong 
repulsion at the typical C2 separation for pairs of Bz's in 
the P configuration. 



C. Friedel lattice (CI) 

The asymptotic form of the pairwise indirect-electronic 
interaction, shown in Fig. [6| predicts a second minimum 
of only 8 meV, comparable to the calculated magnitude 
of surface-corrugation effects. However, in forming a lat- 
tice, there are multiple pairs contributing, effectively in- 
creasing the impact of this effect by a factor of 3. Again, 
assuming a hexagonal lattice and using A = 0.25, we get 
the experimental minimum of a = 10.2 A and a binding 
energy of 3 x 8meV = 24meV. 

The overlayer lattice parameter in phase CI might be 
also influenced by the surface corrugation, since the over- 
layer is essentially commensurate with the 4x4 copper 
unit cell, with a lattice parameter of 10.21 A. The soft- 
ness of the binding curves of Fig. [6] supports this idea. 
The diffusion barriers found in this study, however, sug- 
gests that it is hard to account for the stability of the 
Cl-phase without including the surface- mediated inter- 
actions. Even if the scattering properties of the ben- 
zene adsorbate are closer to those of copper adsorbates 
(A = 0.08, 6 = 0.37r), and therefore does not deter- 
mine the overlayer- lattice parameter, the resulting bar- 
rier should contribute to the stability of this phase. 



V. CONCLUSIONS 

The low barrier to diffusion implies that benzene ad- 
sorbates can respond to weak interactions. The low bar- 
riers could, more generally, allow acene adsorbates to 
self-organize in response to either (or both) direct vdW 
interactions or (and) indirect surface-state-mediated in- 
teractions. There are several conditions for this to occur. 
The adsorption of the acene must be sufficiently weak 
that the process does not destroy the metallic surface 
state yet strong enough that the interaction it mediates 
can account for the regularity of the honeycomb network. 
For benzene molecules we have shown that these condi- 
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tions are present and that the hexagonal-lattice ordering 
found in the CI and C2 phases can be interpreted as 
arising from substrate-mediated and vdW binding, re- 
spectively. 
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